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Abstract 

Mechanically induced unfolding of passive crosslinkers is a fundamental biological phenomenon encountered across 
the scales from individual macro-molecules to cytoskeletal actin networks. In this paper we study a conceptual model 
of athermal load-induced unfolding and use a minimalistic setting allowing one to emphasize the role of long-range 
interactions while maintaining full analytical transparency. Our model can be viewed as a description of a parallel 
bundle of N bistable units confined between two shared rigid backbones that are loaded through a series spring. We 
show that the ground states in this model correspond to synchronized, single phase configurations where all individual 
units are either folded or unfolded. We then study the fine structure of the wiggly energy landscape along the reaction 
coordinate linking the two coherent states and describing the optimal mechanism of cooperative unfolding. Quite 
remarkably, our study shows the fundamental difference in the size and structure of the folding-unfolding energy 
barriers in the hard (fixed displacements) and soft (fixed forces) loading devices which persists in the continuum limit. 
We argue that both, the synchronization and the non-equivalence of the mechanical responses in hard and soft devices, 
have their origin in the dominance of long-range interactions. We then apply our minimal model to skeletal muscles 
where the power-stroke in acto-myosin crossbridges can be interpreted as passive folding. A quantitative analysis of 
the muscle model shows that the relative rigidity of myosin backbone provides the long-range interaction mechanism 
allowing the system to effectively synchronize the power-stroke in individual crossbridges even in the presence of 
thermal fluctuations. In view of the prototypical nature of the proposed model, our general conclusions pertain to a 
variety of other biological systems where elastic interactions are mediated by effective backbones. 


1. Introduction 

In contrast to inert matter, distributed biological systems are characterized by hierarchical network architectures 
with domineering long-range interactions. Even in the absence of metabolic fuel this leads to a highly nontrivial 
cooperative mechanical behavior in both statics and dynamics. Passive collective effects are usually revealed through 
synchronized conformational changes interpreted here as generic folding-unfolding transitions. The experiment shows 
that such systems exhibit coherent macroscopic hopping between folded and unfolded configurations resisting the 
destabilizing effect of finite temperatures (Dietz and Rief, 2008; Thomas and Imafuku, 2012; Erdmann et al., 2013). 

A minimal mechanical model showing the cooperative behavior is a parallel bundle of bistable units linked by two 
shared backbones and its most natural biological prototype is a muscle half-sarcomere undergoing the power-stroke 
(Caruel et ah, 2013). Another example is the unbinding of focal adhesions with individual adhesive elements coupled 
through a common elastic background (Erdmann and Schwarz, 2007). Similar behavior has been also associated with 
mechanical denaturation of RNA and DNA hairpins where the long-range interactions are due to the prevalence of 
stem-loop structures (Liphardt et ah, 2001 ; Bosaeus et ah, 2012). Other molecular systems with cooperative unfolding 
include protein /3-hairpins (Munoz et ah, 1998) and coiled coils (Bornschlogl and Rief, 2006). The backbone domi¬ 
nated internal architecture in all these systems leads to a mean-field type mechanical feedback which is also ubiquitous 
in bi-stable social systems (Kometani and Shimizu, 1975; Desai and Zwanzig, 1978). 


Email address: matthieu. caruel@u-pec. fr (M. Camel) 
Preprint submitted to Journal of the Mechanics and Physics of Solids 


January 8, 2015 





In this paper we systematically study the mechanics of the minimal model in the setting which can be directly asso¬ 
ciated with skeletal muscles. We recall that the mechanically induced conformational change (power-stroke) in skele¬ 
tal muscles takes place in myosin heads (cross-bridges) that are bound in parallel to actin filaments (Smith et ah, 2008; 
Linari et ah, 2010; Guerin et ah, 2011; Erdmann and Schwarz, 2012; Piazzesi et ah, 2014). The (thermo)mechanical 
behavior of this system was first analyzed by Huxley and Simmons (Huxley and Simmons, 1971) who interpreted 
the pre- and post-power-stroke conformations of the myosin heads as discrete chemical states (spin model). Similar 
ideas have been independently advanced in the studies of bistable adhesion clusters (Bell, 1978) and in other applica¬ 
tions ranging from lahn-Teller effect and ripples in graphene sheets (Bonilla et ah, 2012) to unzipping of biological 
macromolecules (Gupta et ah, 2011; Prados et ah, 2012). 

Since in (Huxley and Simmons, 1971) the behavior of the spin model was studied only in a hard device (pre¬ 
scribed displacements), the cooperative effects were not found. To understand this surprising result we study the 
zero-temperature analog of the Huxley and Simmons model. We show that for this system the structures of the energy 
landscape in hard and soft (prescribed forces) loading conditions are rather different. In particular, we explain why 
in this model the collective behavior at finite temperature can be expected in the soft but not in the hard device. To 
capture the coherent hopping in the hard device case, we regularize the spin model in two ways. First, following 
(Marcucci and Truskinovsky, 2010) we replace the discrete chemical states (hard spins) by a continuous double-well 
potential with a finite energy barrier (snap-spring model). Second, to take into account the elastic interactions be¬ 
tween individual crosslinkers, we introduce a series spring mimicking the backbone elasticity (Wakabayashi et ah, 
1994; Huxley et ah, 1994) and bringing in mean-field interactions. We show that in the snap-spring model the energy 
barriers separating the synchronized states are still markedly higher in a soft device than in a hard device which pro¬ 
vides an explanation for the observed retarded relaxation in isotonic experiments on skeletal muscles (Reconditi et ah, 
2004; Piazzesi et ah, 2002; Decostre et ah, 2005). 

An interesting peculiarity of the snap-spring model is that the relaxed potential, representing the global minimum 
of the energy, is always convex in a soft device but is only convex-concave in a hard device. This means that the 
macroscopic stiffness, which is always positive in a soft device, can become negative in a hard device. The negative 
stiffness (metamaterial) response (Nicolaou and Motter, 2012), which persists in the continuum limit, clearly contra¬ 
dicts the intuition developed in the studies of systems with short-range interactions. The non-convexity of the ground 
state energy in a hard device means that the system is non-additive and cannot be relaxed through the mixing of folded 
and unfolded units. 

To illustrate our general results, we consider in some detail the special case of skeletal muscles where we can make 
quantitative estimates by using realistic parameters. Our analysis shows that the height of the microscopic energy 
barriers for the power-stroke in individual cross-bridges is of the order of the energy of thermal fluctuations. This 
implies that the cross-bridges can undergo conformational changes in a non-cooperative stochastic manner. However, 
the presence of long-range interactions creates a bias in the individual folding-unfolding equilibria in the form of a 
macroscopic barrier on top of which the microscopic barriers are superimposed. We call this barrier macroscopic 
because its height is proportional to the number of elements in the system. Due to the presence of the macroscopic 
barrier, the folding transitions in individual cross-bridges become energetically preferable only after the top of this 
barrier has been reached which means that individual cross-bridges have been synchronized. These observations 
suggest that the elementary contractile unit of skeletal muscles has evolved to control the state of a large assembly of 
folding elements through a mean-field type mechanical feedback. Due to such passive synchronization, the power- 
stroke takes place collectively which obviously amplifies the mechanical effect. 

While we focus in this paper on the athermal behavior, our analysis reveals the origin of the anomalous thermo¬ 
dynamics and kinetics of skeletal muscles and similar systems observed at finite temperatures (Caruel et ah, 2013). A 
detailed study of the temperature effects on the collective unfolding will be presented elsewhere. 

The paper is organized as follows. In Section 2, we study the equilibrium mechanical behavior of the spin model 
and show that already in this minimal setting the behavior of the system in soft and hard devices is different. The snap¬ 
spring model is introduced in Section 3 where we demonstrate that it removes the degeneracies of the spin model and 
effectively interpolates between the soft and hard device behaviors. In the same section, we also study the fine structure 
of the energy landscape separating the coherent states of the system and introduce a reaction coordinate to describe 
the successive individual folding-unfolding transitions constituting the collective unfolding. The adaptation of the 
snap-spring model for skeletal muscles is presented in Section 4. Finally in Section 5 we present our conclusions. 
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2. The hard spin (HS) model 

Consider the behavior of an elementary cluster of N bistable units connecting two rigid backbones. In the spin 
model (Huxley and Simmons, 1971), each crosslinker is represented by a bistable potential connected to a series 
spring. The potential, representing two folding configurations, is assumed to have infinitely narrow energy wells 
representing two chemical states, see Fig. 1. The potential describing individual spin units can be written in the form 


«hs (•*') = 


jvo 

1° 


if x = 0, 
if x = —a. 


( 1 ) 


Here the spin variable x takes two values, 0 and - a , describing the unfolded and the folded conformations, respectively. 
By a we denoted the “reference” size of the conformational change interpreted as the distance between the two energy 
wells. With the unfolded state we associate an energy level vo while the folded configuration is considered as a zero 
energy state. The potential (1) is shown schematically in Fig. 1(a). 
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Figure 1: Hard spin model of a parallel bundle of bistable crosslinkers, (a) Energy landscape of an individual crosslinker; (b) N crosslinkers loaded 
in a soft device. 


In addition to a spin unit with energy (1) each cross-bridge contains a linear shear spring with stiffness k\ see 
Fig. 1(b). The energy of the elastic spring is u E (x) = kx 2 /2 and the energy of the whole crosslinker is 


u = u HS (x) + u E (z - x). 


( 2 ) 


Without loss of generality, we can assume that the reference length of the linear spring is already incorporated into 
the definition of the elongation z. Notice that the mechanical system with energy (2) has a multi-stable response, see 
Fig. 2. 

To non-dimensionalize the resulting model, which one can associate with the names of Huxley and Simmons 
even though they never considered such parallel connection explicitly (for this representation, see Marcucci and 
Truskinovsky (2010)), we choose a as a characteristic length, associate the characteristic energy scale with ka 2 and 
normalize forces by ka. The only remaining dimensionless parameters of the model are N and vo and we can write 
the dimensionless energy of the system (per crosslinker) in the form 



(3) 


Here, for convenience, we preserved the same notations for non-dimensional quantities. 

In a hard device each crosslinker is exposed to the same total elongation z and thus the individual units are 
independent. In the soft device case, where the control parameter is the total tension T, the energy per crosslinker is 


w(x,z;t) = v(x,z) -tz = 


i N 

iy 

N 


(1 +xi)v 0 + -(z-XiY -tz 


(4) 
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Figure 2: Behavior of a single crosslinker in the spin model, (a) Equilibrium states for various z\ (b) Corresponding tension levels. Dashed lines, 
metastable states; bold line, global minimum. Arrows show the response to sudden shortening including a frozen elastic phase (A —» B i) and a 
subsequent phase equilibration ( B\ — > B 2 ). Here we used vq = 0. 


where t — T/N is the force per crosslinker. Now for each crosslinker both x, and z are internal degrees of freedom and 
the individual units can no longer be considered as independent. Indeed, if we minimize out the global continuous 
variable z by solving dw/dz \ t _(*.} = 0, we obtain 


" t+ N 


1 N 


(5) 


which after inserting in Eq. 4 shows that the partially minimized energy depends quadratically on Yj'iLi x i > 




x? 


— - tXj + (1 + X;)v 0 - — 




One can see that the transition from hard to soft device introduces a mean-field interaction among the crosslinkers 
which, as we show later, is ultimately responsible for the cooperative behavior. 


2.1. Mechanical equilibrium in a hard device 

To describe the equilibrium response of the HS model in a hard device, we need to compute the local minima of 
the mechanical energy (3) at fixed z. Since each of the internal degrees of freedom jc,-, for 1 < i < N can take only 
two discrete values, x,- = xo = 0 and x, = x\ = —1, a given equilibrium state is characterized by the distribution of 
the N crosslinkers between the two spin configurations. Due to the permutational invariance of the problem, each 
equilibrium state is fully characterized by a discrete order parameter representing the fraction of crosslinkers in the 
folded state, 

1 N 

1=1 

where a, — 1 if x, — — 1 and a,- = 0 if x,- = 0. 

In Appendix A we show that all equilibrium configurations of this type correspond to local minima of the energy 
(3). At a given value of the order parameter p, the energies of all such metastable states are equal to 

v(p; Z) = p ^ (z + l) 2 + (1 - p)(^Z 2 + v 0 ). (7) 

This is a linear combination of the energies of two limiting configurations, one fully folded with p = 1 and the energy 
|(z + l) 2 and the other one fully unfolded with p — 0 and the energy \z 2 + vo. The absence of a mixing energy is a 
manifestation of the fact that the two coexisting populations of crosslinkers do not interact. 
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Figure 3: Mechanical response of the HS model in a hard device with N = 5. (a) Energy levels of the metastable states (p = 0, A,..., 1) for 
different applied elongations, (b) corresponding tension-elongation relations, (c) and (d) are details of (a) and (b) with an illustrated response path 
to a fast loading experiment including a frozen elastic phase (A —> Bi) followed by a subsequent phase equilibration ( B\ —» B 2 ). Thick lines, global 
minimum corresponding to p = 0 (resp. p = 1) for z > z* (resp. z < z*) • f* = vq and z* = vq - 1/2. Parameters are vq = 1 and N = 5. 


The tension-elongation relations along metastable branches parameterized by p can be written as 

Q 

Kp\ z) = — v(p; z) = z + p, 
oz 


( 8 ) 


so for given p we obtain equidistant parallel lines, see Fig. 3. A peculiar feature of the HS model is that the domain 
of hysteretic behavior extends indefinitely because the spin system does not have any stress thresholds. 

To find the global minimum of the energy we need to perform at each value of z an additional minimization over 
the discrete variable p. If we compute the derivative 


d_ 

dp 


V(P\ Z) = Z + - - V 0 . 


which does not depend on p, we obtain that for z > Z„ where z, = v ; o - 1/2, the global minimizer is the folded state 
with p = 1 and for z < z» it is the unfolded state with p = 0, see Fig. 3. The global minimum energy profile exhibits a 
kink near the crossing (folding) point. Similar kinks associated with unfolding of hairpins and other folding patterns 
have been observed in the energy profiles reconstructed from single molecule force spectroscopy measurements of 
proteins and nucleic acids (Gupta et ah, 201 1). 

Observe that the ground state energy in this model is nonconvex independently of the number of units. This means 
that the energy is not convexified through the formation of mixtures as in systems with short-range interactions (see, 
for instance, Puglisi and Truskinovsky (2000)). The reason is that this mean-field system is strongly non-additive and 
all mixed states are energetically unfavorable due to high cost of mixing. Somewhat similar situation takes place in 
theory of elastic phase transitions where the relaxed energy is also nonconvex in the general case (it is only quasi- 
convex) which is again the consequence of long-range elastic interactions exemplified by the gradient constraint (see, 
for instance. Ball (2002)). 

A consequence of the energy nonconvexity is the non-monotonicity of the force-elongation relation shown in 
Fig. 3(b,d). More specifically, the system exhibits a negative stiffness at the point where all crosslinkers collectively 
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flip from the folded to the unfolded state. Similar mechanical behavior has been recently artificially engineered in 
metamaterials by drawing on the Braess paradox for decentralized globally connected networks (Cohen and Horowitz, 
1991; Nicolaou and Motter, 2012). As in our case, the mean-field type coupling in metamaterials is achieved via 
parallel connections with multiple shared links. Biological examples of systems with negative stiffness are provided 
by unzipping RNA and DNA hairpins (Woodside et ah, 2008; Bosaeus et ah, 2012). 


2.2. Mechanical equilibrium in a soft device 

Consider now the HS model loaded in a soft device when the equilibrium states correspond to local minima of the 
mechanical energy (4). An equilibrium state is again fully characterized by the fraction of units in the folded state p, 
defined by Eq. 6. We can then write the (marginal) energy of the partially equilibrated system as 


w(p,z\ t ) = v(p,z) - tz. 


where v is the energy of the metastable branch parametrized by p in the hard device case, see Eq. 7. 

Then we eliminate z by using Eq. 5. Each value of the order parameter p defines a branch of local minimizers of 
the energy (4) parameterized by t, see Appendix B. At a given value of p, the energy of a metastable state reads 


w(p; t) 


1 , 1 
--r + pt+ -p( 1 


■p) + (1 - P)v o- 


(9) 






Figure 4: Mechanical response of the HS model in a soft device, (a) Energy levels of the metastable states (p = 0, j,..., 1) for different applied 
forces, (b) Corresponding tension-elongation relations. In (c) and (d) we zoom into domain near t = t* = vq and show schematically the response 
of the system to a sudden application of a load increment with an elastic phase (A —» C\) followed by a folding phase C\ —» C 2 . Thick lines, global 
minimum corresponding to p = 0 (p = 1) for t > t* (t < t *). Parameters are, vq = 1 and N = 5. 


In contrast to the case of a hard device, here there is a nontrivial coupling term p( 1 - p) describing the energy of 
mixing. The presence of this term is a signature of a mean-field interaction among individual crosslinkers. Indeed, if 
one element changes configuration, its contribution to the common tension changes accordingly and the other elements 
must adjust to maintain the force balance. The tension-elongation relation associated with a set of metastable states 
sharing the same value of the parameter p can be written in the form 


zip', t) = - 


d_ 

dt 


w(p', t ) 


= t-p. 
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At a given value of p this relation is identical with its counterpart in the case of a hard device, see Eq. 8. The globally 
stable states can be found by minimizing (9) over p. Since d 2 w/dp 2 = — 1, this function is concave in p. Therefore, 
the global minimum is again attained either at p - 1 or p = 0. The energies of the two coherent configurations with 
p = 1 and p — 0 coincide when t - t, - i’ () . 

Our Fig. 4 illustrates the structure of the energy-tension and the tension-elongation relations corresponding to 
different values of p for the system with N — 5. Notice that in contrast to the case of a hard device, the force- 
elongation relation characterizing the global minimum in a soft device exhibits a plateau indicating a discontinuity in 
elongation as the crosslinkers switch collectively at t = t„ from unfolded to folded conformation. 

Notice that the plateau replaces the region of negative stiffness detected in the hard device case and the force- 
elongation relation becomes monotone. This shows that even in the continuum limit the stable “material” responses 
of our system in hard and soft devices remain different. Such behavior would be completely unexpected from the 
perspective of classical statistical mechanics, however, it is characteristic of systems with domineering long-range 
interactions (Dauxois et al., 2003). 

2.3. Energy landscape 

As we have seen in the previous sections, the globally stable state of the HS system correspond to one of the 
two coherent configurations characterized by p = 1 and by p = 0. We can now pose the question about the size of 
the energy barrier separating these two configurations. To access the energy barriers and to find the transition states 
(saddle points of the energy) we study the energy dependence on p. For general values of the loading parameters this 
dependence was found to be linear in a hard device, indicating that there is no conventional barrier, and concave in a 
soft device which means that there is a potential energy barrier. This observation shows that a switch from unfolded 
to folded configuration in a soft device carries an energetic cost while in a hard device the transition is cost-free. 


C + C B + B A + A 




Figure 5: (a) Energy landscape at the global minimum transition for the HS model. Solid lines, hard device at z = Z*\ Dashed lines, soft device at 
t = f*. Dots represent the energy of the different configurations for a system with N = 5; lines correspond to the limit N — * oo. Parameters are 
vo = 1 and N = 5. (b) Representation of the behavior of a system with two crosslinkers with vo = 0 imposing ?* = 0 and z* = -1/2, the transition 
of the global minimum tension-elongation curve (thick line) occurring in a stress-free configuration in a soft device. Dashed lines, metastable states 
p — 0 and p = 1. The intermediate stress-free configuration is obtained either by mixing the two geometrically compatible states B and D in a hard 
device which results in a B + D structure without additional internal stress or by mixing the two geometrically incompatible states A and C in a soft 
device which results in a A + C structure with internal residual stress. 


The (collective) transition takes place at z = Z* in a hard device and t — t* in a soft device. The difference between 
the corresponding energy landscapes at the threshold values of parameters is illustrated in Fig. 5 (a). For the ease of 
comparison the energy minima are shifted to zero in both loading conditions. Each black dot represents the energy 
level of a particular metastable state and the dotted line represents the set of such metastable states in the continuum 
limit N —» oo when p becomes a continuous variable. The barriers separating individual metastable states are not 
defined in the spin model. A simple way to recover microscopic barriers by switching from hard to soft spins is 
discussed in Sec. 3; for another regularization approach see (Benichou and Givli, 2013). 
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To understand the origin of this peculiar behavior of the energy landscape it is instructive to consider the minimal 
HS system with N - 2\ see Fig. 5 (b). Here for simplicity we assumed that i'o = 0 implying 1, = 0 and z* = -1/2. 
The two “pure” configurations are labeled as A (p = 0) and C (p — 1) at f = f. = 0 and as D (p = 0) and B (p — 1) at 
z — z* = -1 /2. In a hard device, where the two elements do not interact, the transition from state D to state B at a given 
z = z* goes through the configuration B + D which has the same energy as configurations D and B: the crosslinkers 
in folded and unfolded states are geometrically perfectly compatible and their mixing requires no additional energy. 
Instead, in a soft device, where individual elements interact, a transition from state A to state C taking place at a given 
t = 0 requires passing through the transition state A + C which has a nonzero residual stress. Individual crosslinkers 
in this mixture state have different values of z and therefore the energy of the stressed “mixed” configuration A + C is 
larger than the energies of the “pure” unstressed states A and C. We conclude that the macroscopic barrier in a soft 
device is higher than in a hard device because in a soft device a transition is a genuinely cooperative effect requiring 
essential interaction of individual elements while in a hard device the conformational change in different units takes 
place independently. 



Figure 6: Summary of the behavior of the HS model. The main frame shows the tension-elongation relations corresponding to p = 0,0.1,..., 1 
for a system with N = 10 (gray lines) and the tension-elongation relation in the global minimum (thick lines). The horizontal dotted lines show 
the limit of the domain where a macroscopic energy barrier is present in a soft device. The satellite frames (above, hard device; right, soft device) 
show the energy barriers corresponding to various transitions (A—>B, C—>D,...) shown by dashed lines in the main frame. 


For general values of the loading parameters, we consider the N —> °o limit where the energy landscape becomes 
smooth and the barrier is located at the saddle point p = /?„, where p t = t — v o + 1/2 is the solution of dw/dp = 0. 
Since 0 < p* < 1, we obtain that the macroscopic barrier and the ensuing cooperative effects exist for < t < t+ with 
t- = t, ~ 1/2 and t+ = t* + 1/2, in the soft device setting. Notice that the boundaries of this interval correspond to the 
threshold f(0;z„) and t( 1; ) delimiting the region with negative stiffness in the hard device case, see Fig. 3. 

We summarize the results obtained within the HS model in Fig. 6. In the main frame, we show the tension- 
elongation relations corresponding to p = 0,0.1,.... 1, for a system with ;V = 10 (gray lines). Thick lines show the 
tension-elongation relations in the global minimum characterized by a transition located at z* in a hard device and t, 
in a soft device. The horizontal dotted lines show the limits of the domain where a macroscopic barrier is present in a 
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soft device. On the satellite frames (top - for a hard device and right - for a soft device) we show the energy barriers 
between the homogenous configurations in the case N —» oo (solid lines) and for discrete values of p (dots) at different 
values of the loading. 

To conclude, we have shown that the HS model exhibits different mechanical responses in a hard and a soft 
devices even though in both cases the global minimum of the energy is achieved for homogenous configurations 
This behavior originates from the presence of long-range interactions between the crosslinkers. These interactions 
introduce an elastic feedback in the soft device case which creates a macroscopic energy barrier for the cooperative 
folding-unfolding process. 


3. The snap-spring model 

The minimal description of unfolding transition provided by the spin model has its limitations. For instance, the 
description of the barriers in the HS model is incomplete because it does not take into account the microscopic energy 
barriers between the states with different values of the order parameter p. Another problem is that in the hard device 
setting the spin model is degenerate because configurations with different values of the order parameter are equivalent. 
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Figure 7: Snap-spring model of a cluster, (a) Dimensional energy landscape of a bistable crosslinker, (b) Structure of a parallel bundle containing 
N crosslinkers in a soft device. 


To deal with these problems we regularize the spin model by introducing two additional physical mechanisms. 
First, following (Marcucci and Truskinovsky, 2010) we replace hard spins by snap-springs also known as soft spins, 
so that x becomes a continuous variable. For simplicity we assume that the corresponding double-well potential can 
be represented as a minimum of two parabolas, see Fig. 7(a). By using non-dimensional variables we can then write 

. . _ J %(x) 2 + v 0 if x > / 

UssW - |i fei(jc+ 1} 2 ifx < z 

Here 1 is the dimensionless position of the energy barrier, v ; o is the dimensionless energy bias of the unfolded state and 
k\ and &q, are dimensionless elastic moduli of the folded and unfolded states, respectively. Interestingly, a comparison 
with the reconstructed potentials for unfolding biological macro-molecules shows that this approximation may be in 
fact very good (Gupta et al., 2011). A spinodal region can be obviously added to the potential of the bi-stable unit, 
however, in this case we lose transparency without gaining essential effects. 

The important observation is that while in the snap-spring model the bottoms of the energy wells remain the 
same as in the spin model, the barrier separating the two conformational states is now well defined, see Fig. 7. The 
mechanical response of a single crosslinker with an attached series spring is governed by the dimensionless energy 

1 

u = u ss (x) + -(y - x) , 

where y is the total elongation. The solutions of the equilibrium equation u' ss (x) = y-x are shown in Fig 8. Notice the 
multi-valuedness of the relation linking the variables x and y and the appearance of the spinodal branch x* connecting 
the two stable branches xi(y) and xo(y). 
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Figure 8: Behavior of a single crosslinker in the snap-spring model, (a) Equilibrium positions for various y; (b) Corresponding tension levels. 
Solid lines, metastable states; dashed lines, unstable state; bold line, global minimum. Arrows indicate the response to a sudden shortening with an 
elastic unloading in the unfolded state (A —» B\) followed by the conformational change to the folded state ( B\ —» B 2 ). Parameters are, A\ = 0.4, 
A 0 = 0.7, / = -0.3 


The mechanical independence of the crosslinkers in a hard device disappears if we take into account the finite 
stiffness of backbone which, in the case of skeletal muscle, corresponds to the combined stiffness of actin and myosin 
filaments (Wakabayashi et ah, 1994; Huxley et ah, 1994; Ford et ah, 1981; Mijailovich et ah, 1996; de Gennes, 2001). 
Following (Jiilicher and Prost, 1995), we use a lump description of backbone elasticity by introducing an additional 
elastic spring with stiffness Ab = kb/(Nk) and the energy Ub(x) = NAbX 2 / 2. If we attach this spring in series to our 
parallel bundle of crosslinkers, see Fig. 7(b), we can write the total energy of the system per crosslinker in the form 


1 

v(*,v;z) = -2 

(=1 


Us S (xd + ^(y - Xij 1 + y (z - yj 1 


( 10 ) 


In the hard device case z is the control parameter, x, are the continuous microscopic internal variables generalizing 

the spin variables in the HS model and y is a new continuous mesoscopic internal variable. Notice that now even in 

a hard device the individual crosslinkers are not independent; the implicit mean-field interaction becomes obvious if 

the variable y is adiabatically eliminated (minimized out) by solving =0. We obtain 

°y i z.x 


y(x; z) 


1 


1 N A 

AbZ+ N^ Xi 

i= 1 / 


1 + Ab 

which after inserting into Eq. 10 gives the partially equilibrated energy 
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Observe that the quadratic term in £ x,- vanishes when the elasticity of the backbone becomes infinite (Ab —> oo) 
showing that in this limit the long-range interactions disappear. 

The resulting snap-spring model can be viewed as a regularization of the HS model. To recover the HS model in 
a hard device case from Eq. 10, we need to perform the double limit: k t n —» °° and Ab —> oo. The first of these limits 
ensures that x becomes a spin variable while the second guarantees that y = z- To obtain the HS model in a soft device 
we need to consider the triple asymptotics: k\ o —> oo, Ab —> 0 and z —> °o where the last two limits must be linked in 
the sense that AbZ —> t which ensures that the force per crosslinker t remains finite. 

In a soft device, the total energy per crosslinker in the snap-spring model can be written as 

Uss(xi) + ^(y - X,) 2 + y (z - yf - tz , 

where t = TIN is again the applied force per crosslinker. 


w(x,y,z\ t) = v(x,y,z) - tz = ^ ^ 
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3.1. Mechanical equilibrium in a hard device 

To find equilibria in a hard device we need to solve the following system of equations 


dv 


dxj 


= 0 for all 1 < i < N 


dv 


dy 


= 0 


z,M 


Equation 12a have up to 3 solutions that can be parameterized by y. 


(12a) 

(12b) 


x\ iy) = (1 - di)y- Au iix t < l 

'X 0 (y) = (l-Ao)y, if Xj>l 

x* = l 


(13) 


where we redefined the dimensionless parameters 


d 0 = 


_k o_ 

l+^o 


A, 


h 

\ + k\ 


The solution of Eq. 12b given by Eq. 11 allows to express the mesoscopic variable y through the microscopic variables 
X/. Because of the permutational invariance, the equilibrium solution of Eq. 12 is fully characterized by the fraction p 
of crosslinkers in the folded conformation xi (v) and the fraction r of crosslinkers in the unfolded conformation To (v)- 
The fraction of crosslinkers in the “spinodal point” x, is then q — 1 - p — r. 

Using (11), we can eliminate the variable y and obtain an explicit representation of the microconfiguration in terms 
of ( p,q,r ): 


y{p,q,r\z ) = 
X| (p, q, r, Z) = 


A-b 


- 

,q,r)\ 


A b + A (p 

1-di 

A b + A(p, q, r) 


pA\ — ql \ 
A b I 


(AbZ-pAi + ql ) - Au 


x 0 ip, q, r, z) = - 7 —-- (AbZ - pA\ + ql) . 

A b + A (p, q, r) 


(14) 

(15) 

(16) 


Here A( p, q,r) = pA \ + q + rA () represents the equivalent stiffness of the parallel bundle of crosslinkers in a mixed con¬ 
figuration parameterized by (p, q, r). The energies of the equilibrium configurations can be now computed explicitly. 
For a given ip, q, r) we obtain 


vip,q,r,z ) 


A b \pA\iz + l) 2 + qiz - l) 2 + rA Q z 2 ] + pA\ [rd 0 + qi 1 + 2/)] - ql 2 iq + A b ) q f 


2 iA b + A ip,q,r)) 
The corresponding tension-elongation curves can be written as 


2(1-do) 


■iq + r) I'o. (17) 


tip, q, r, z) 


—vip,q,r\z) 

dz 


A b Kjp, q, r) I pAi - ql \ 
A b + A(p,q,r) \ A ip,q,r)j 


Each triple ip,q,r) defines an equilibrium branch which extends between the two limits [z m r, z sup ] induced by the 
inequalities X| < / and xo > /. We obtain 


Z m tip, q, r) 


Zsupip, q, r) 


l [A b + A (p, q, r)] + (1 - A Q )ipA\ - ql) 

Abi 1 - d 0 ) 

U + di) [A b + A jp, q,r)\ + j 1 - Ai)jpA x - ql) 
di(l-d!) 


( 18 ) 
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Figure 9: Mechanical response of the snap-spring model in a hard device, (a) Energy levels of all the (p , q, r ) configurations for the case N = 3 at 
different applied elongations, (b) corresponding tension-elongation relations. Solid lines, metastable states with p — 0, 1 /3,2/3, 1 and r = 1 - p, 
q = 0; dotted lines, unstable states with q ± 0; thick lines, global minimum corresponding to p = 0, r = 1, for z > Z*, and to p = 1, r = 0, for/ 
Z < z*. (c,d) blow-up of (a,b) illustrating the response of the system to abrupt shortening with an elastic unloading (A —> B\) followed by a massive 
conformational change in isometric conditions ( B\ —> B 2 ). Parameters are, A\ = 0.4, Aq = 0.7, / = -0.3, Ab = \, and N = 3. 


The analysis presented in Appendix A shows that all equilibria with q = 0 are stable while all the ones with q ± 0 are 
unstable. Therefore, as in the HS model, the metastable configurations in the snap-spring model can be parameterized 
by a single parameter p. The obtained results are illustrated in Fig. 9. 

We now show that the global minimum of the energy is again achieved exclusively on homogeneous configurations 
(1,0,0) and (0,0,1). Assuming that q — 0, r — l-p and computing the second derivative of v in (17) while interpreting 
p as a continuous variable leads to 


<9 2 1 [di/li (z + 1) + do (Ti - 

- —zV(p\Z) = -o- 

d P~ 2 [pAi + (l - p) A 0 + A b y 


(19) 


This inequality shows that the energy is concave, which means that the ground states are necessarily synchronized 
and separated by mixed configurations with higher energy levels forming a macroscopic energy barrier. Observe that 

o2 

for Ab —> oo we have J^v(/?;z) —> 0 and thus the macroscopic barrier disappears and we recover the degeneracy of 
the HS model. The switch between the two homogeneous states takes place at the elongation z = z» which solves 
v(0,0, l;z„) = v(l,0,0;z«). We obtain 


Z.* 


( [(To - Ti) Ab] 1 [Ti (Tq + T/,) - (To + T*) (Ti + T/,) | 
< Ti + Ab 1 
i TjT* l '° ~ 2 


if Ti * T 0 
if T, = T 0 


where p = ToTi + 2(Ti - To)v'o > 0. Notice that we recover the position of the transition point of the HS model, 
z, = Vo - 1/2, when considering the symmetric case with Ti = 1 and Ab —> oo. 

While the global minimum path in the snap-spring model has the same structure as in the HS model, we see that at 
the transition point z = z* the energies of the mixture states are now strictly higher than the energies of the coexisting 
pure states, see Fig. 9(a,c). 
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The ensuing force-elongation relations, presented in Fig. 9(b,d), show that the singular metamaterial behavior ex¬ 
hibited by the HS model in a hard device is not regularized in the snap-spring model where the stiffness corresponding 
to the globally stable response is still equal to minus infinity at the transition point. 


3.2. Mechanical equilibrium in a soft device 

To find equilibrium states in the snap-spring model loaded in a soft device we need to solve the system 

dw 

= 0 for all 1 < i < N 

= 0 

= 0 


dx. 



\ x i*i 

dw 


~dy 


dw 


dz 



As in the hard device case, each crosslinker can be in three states and the equilibrium branches can be parameterized 
by the triplet (p,q,r). After elimination of the internal degrees of freedom x, and y, the corresponding partially 
equilibrated energy can be written in the form 

w(p, q, r, z; t) = v(p, q, r, z) - tz 

where the function v is given by Eq. 17. After elimination of z, the solution of the full mechanical equilibrium is 
obtained, 


yip, q, n t ) = 


1 ^ pA\ - ql 

A (p, q, r) A (p, q, r) ’ 

xoip, q, n t) = —^r(t - pA\ + ql), 

A ip, q, r) 

xi(p,q,r;t ) = - pAi +ql)-A l . 

A ip, q, r ) 

The energy of a configuration ( p, q, r ) can be again computed explicitly. 


( 20 ) 

( 21 ) 

( 22 ) 


w(p, q, r;t) - — - ( -—t- 


1 


2 \A b A (p,q,r) 


, pA\ - ql pA\rAo - q 2 1 1 + 2pA\ql 
t + - - —t + - 


A (p, q, r) 


2A (p,q,r) 


+ (q + r)v o + 


ql 2 


1 - To’ 


and the corresponding tension-elongation relations read 

d 


z (p, q, r; t) = ~—w{p, q,r\t) = \— + 


dt 


1 


1 


A b h(p,q,r)J A (p,q,r) 


pA\ - ql 


Similarly, we can obtain the lower and upper limits for a branch labeled by (p, q, r), 

tsupip, q, r) = \ + A (p, q, r) + pA\ - ql. 


1 -T, 


(23) 


t m ((p,q,r) = 


1 - To 


A (p, q, r ) + pAi - ql. 


The stability analysis presented in Appendix B shows again that configurations with q ^ 0 are unstable and 
therefore all metastable states can be parameterized by a single parameter p. Since 


d 2 „ 

^ ( P ;r) = - 


A 


(Ti - To) - (t - pAi) + 

\pA\ + (\ - p) T 0 ] 3 pAi+(l-p)A 0 


< 0 , 


the global minimum is again attained either at p = 1 or p — 0 and these pure states are robust and separated by a 
macroscopic energy barrier. Notice that in contrast to the hard device case, see Eq. 19, the presence of this barrier 
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does not depend on the backbone stiffness. The switch between them, signaling a collective folding-unfolding, takes 
place at t — t, where vi>(l,0,0,f*) = vv(0,0,1,?»). We obtain 

U = Wo - T|) 1 [do/li - Vp/Mi] 

which simplifies into 1, = vo if 4 | = Aq . Note that in this later case, the value of t, is the same as in the HS model. 



-1 0 1 



-1 0 


t/t. 


6z = z- z. 


Figure 10: Mechanical response of the snap-spring model in a soft device, (a) Energy levels of all the ( p,q , r) configurations for the case N = 3 
at different applied tensions, (b) corresponding tension-elongation relations. Solid lines, metastable states with p = 0,1/3,2/3,1 and r = 1 - p, 
q = 0; dotted lines, unstable states with q ± 0; thick lines, global minimum corresponding to p = 0, r = 1 for t > U and to p = 1, r = 0, for 
t < t*. (c,d) zoom in of (a,b) illustrating the response of the system to abrupt shortening with an elastic unloading (A —» B i) followed by a massive 
conformational change in isometric conditions (Ci —» C 2 ). Parameters are as in Fig. 9. 


The equilibrium behavior of the snap-spring model in a soft device is illustrated in Fig. 10. As in HS model, the 
globally stable response contains an extended plateau at t = f* where the systems collectively switches between fully 
folded and fully unfolded configurations. The constitutive behavior remains quantitatively the same in the continuum 
limit N —» oo. 


3.3. Energy landscape 

We now study the size of the energy barrier separating the homogenous configurations characterized by p = 0 
and p - 1. In the HS model, we have demonstrated that such macroscopic barrier exists only in the soft device case. 
In the snap-spring model the mechanical feedback introduced by the backbone elasticity introduces a variable degree 
of cooperativity between the crosslinkers which results in an increased energy of the mixed states already in a hard 
device as shown in Fig. 9. 

To find the energy minimizing “reaction path” connecting the homogeneous states through the set of inhomoge¬ 
neous metastable states with 0 < p < 1. The energy of these intermediate state can be made explicit by putting 
q — 0, r = 1-pin (17). We obtain. 


v(p\ Z ) 


1 

2 (/li + A (p)) 


\A h (pA i(z + l) 2 + (1 - p)A 0 z 2 ) + pTi(l - p)/l 0 ] + (1 - p)v o, 


(24) 


with 

Mp) = pA] + (1 - p)A 0 . 
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We have shown that for A/, < °° the function v(p: z) is concave in p, see Eq. 19, which signals the presence of a 
macroscopic energy barrier. However, since the variable p is discrete, this information is incomplete and we need to 
also account for the microscopic barriers separating configurations with different values of p. Such barriers, associated 
with the conformational changes in individual crosslinkers, were essentially infinite in the HS model. 

To reconstruct the fine structure of the energy barriers, we consider a configuration with N\ crosslinkers in the 
folded state, No crosslinkers in the unfolded state and N, = N — N\ -No crosslinkers switching collectively from the 
unfolded to the folded state. The initial stable configuration is fully characterized by the parameter p = N\/N and we 
denote by a — N*/N the fraction of switching crosslinkers which satisfies 0 < a < 1 — p. 

To find the barrier which the system has to overcome, we need to choose a microscopic “reaction path” separating 
the initial configuration characterized by p and the final configuration characterized by p + a. Assuming for simplicity 
that the switching crosslinkers are characterized by the same strain variable x, it is natural to choose x as a “reaction 
coordinate”. 

Due to permutational invariance, the choice of the N t switching crosslinkers is arbitrary, and for commodity, we 
select x\,...,xn,- The energy landscape (per crosslinker) along the chosen reaction path can be recovered if we 
minimize out the rest of the internal variables x,v-,v,,... , x N ,y. We obtain 


v(p, a, x;z) = - 


A b [pA\{z + l) 2 + (1 - p - a)A 0 z 2 ) + pA\(\ - p - a)A 0 


Ab + A (p) + (1 - Ag)a 


u ss (x) + 


■ (1 - P ~ a)v 0 


1 Ab(z - x) 2 + pA \(x + l) 2 + (1 - p - a)At)X 2 


Ab + A (p) + (1 - To)!? 


(26) 


where the first term does not depend on x; it converges to v(p\ z) given by Eq. 24 in the thermodynamic limit when 
N —» oo. The second term vanishes in the thermodynamic limit because it is proportional to a which goes to zero 
when N —> oo. 

For convenience, we map the reaction coordinate to the interval [p, p + a] by replacing x with a stretched variable 
£ defined by 


£ = p + a 


x - x 0 (p; z) 
x\(p + a;z)~ x 0 (p;z)’ 


(27) 


where xq and x\ are the locations of the bottoms of the energy wells defined by Eqs.15 and 16. One can see that the 
variable ( can be viewed as a global reaction coordinate that encompasses the local reaction coordinates of individual 
transitions in a step by step manner and thereby defining the position of the system on the whole reaction path (see 
also Truskinovsky and Vainchtein (2003)). 

The values £ = p and £ = p + a are associated with the metastable states v(p, q — 0,r = 1 - p;z) and v(p + a,q — 
0, r — 1 — p — a; z), respectively, see Eq. 7. At % = p + a(l - Xq)/(x\ - xo), we have x = / and the energy has a local 
maximum, namely v(p, q — a, 1 - p - a). These states belong to the unstable equilibrium branches characterized by 
q > 0 and shown in Fig. 9 by dotted lines. 

The local “microscopic” energy barriers surrounding a given metastable state characterized by a particular p may 
be either in the direction of additional folding (7T_>) or additional unfolding (#<_) of a fraction a of crosslinkers. The 
height of these barriers can be expressed analytically (from Eqs. 24 and 26 ) as follows 


B_,(p, a;z) = N [v{p, a, /; z) - v(p\ z)], 
B^(p, a;z) = N \v(p - a, a, /; z) - v(p\ z)]. 


(28) 


In the thermodynamic limit, we have a = N*/N 0 and the overall microscopic energy barriers take the form 

B^(p:z) ~ a(p\z) - ^b(p;z), 

BUp:z ) ~ a'(p;z) - ^b'(p\z). 


(29) 


The explicit expressions for the coefficients a, b, a' and b' are given in Appendix C. We observe that in the thermody¬ 
namic limit, the height of the energy barriers between different metastable states has a finite limit while the height of a 
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generic barrier per crosslinker (B^/N or /i,//Vj vanishes. Similar results have been previously obtained for chains of 
bistable elements connected in series (Puglisi and Truskinovsky, 2000; Benichou and Givli, 2011; Manca et al., 2013; 
Tshiprut and Urbakh, 2009). 



Figure 11: Energy landscape at the global minimum transition for the snap-spring model with N = 5. (a) hard device at z = z* \ (b) soft device at 
t = u. Solid lines, successive barriers obtained from Eqs. 26 (a) and 31 (b); Dashed lines, continuum limit, N —* oo. Energy minima are arbitrarily 
set to 0 for comparison. Here, Ab = 1 and other parameters are as in Fig. 9. In particular, we have /to =£ A\ which explains why the curve is not 
symmetric. 


To facilitate comparison with the HS model we first limit our attention to the barriers associated with the transi¬ 
tion between the two globally stable coherent states taking place at z = z* and consider only the energy minimizing 
reaction path which is caracterized by successive single crosslinker conformational change (a = 1/N). The resulting 
energy landscape r(£) at z-Z* is shown in Fig. 1 1(a) where we compare two cases, N — 5 (solid line and metastable 
states marked with dots) and N —» oo (dashed line). At finite N we see the macroscopic barrier , not captured by the 
HS model and a superimposed set of microscopic barriers representing “lattice pinning”. These microscopic barriers 
are due to the discreteness of the problem and they disappear in the continuum/thermodynamic limit N = oo when 
(r —> P- 


In the case of a soft device, the energy of metastable states is obtained by inserting q — 0 in Eq. 9. We obtain 

1 


w(p: t) = \ — + „— 

2 A (p) 


pAx pAiil - p)A Q 

- - 1 + -x- + (1 - P)\’Q. 


A ip) 


2Mp) 


(30) 


The partially equilibrated energy at fixed t, p, a and x is then given by 


1/1 

w(p,a,x\t ) = — - -— I- 


1 


2 \J-b A(p) + (1 - A{))a 


t 2 + P' 0(1 ~ P ~ g)T 0 + 2pA\t + 


2(A (p) + (1 - A 0 )a) 
+ a 


p - a)v o 

+ (1 - p — a)Aox 2 — 2tx 


1 pAi(x + l) 2 

Uss(X) + - -^- 

2 K(p) + (1 - A 0 )a 


(31) 


We again map the reaction coordinate x in the interval [p,p + a] using Eq. 27 where xo(p~,z ) and x\(p\z) have to be 
replaced by Jco {p\ t) and Jcj (p; t ) given by Eq. 21 and Eq. 22, respectively. The microscopic energy barriers for folding 
(B_,) or unfolding ( B, _) of A* crosslinkers are obtained analytically from Eqs. 30 and 31, 

B^(p, a\t) = N [w(p, a, x; t ) - w(p; f)], 

B^(p, a;t) — N \w{p - a, a, x\ t ) - w(p; f)] ■ 

As in the hard device case, we obtain in the thermodynamic limit. 


(32) 


B^(p\ t) ~ dp: t) - jjdip: t ), 

BUp: t) ~ c'(p: t) - jjd'iP* t). 


( 33 ) 
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Figure 12: Energy landscape at the global minimum transition for the snap-spring model in a hard device at different values of the coupling 
parameter A/, with N = 20. Bold line, At, —> oo; solid line, Ab = 1; dotted line, soft device limit as the hard device case where At, —» 0, z* —> oo. 
Energy minima are arbitrarily set to 0 for comparison. Other parameters are as in Fig. 9. 


The explicit formulas for the coefficients c, c', d and d' are given in Appendix D; as one can expect, these coefficients 
do not depend on the backbone stiffness A b . Notice that the result is fully analogous to what we obtained in the case 
of a hard device, see Eq. 29: in the thermodynamic limit, the individual overall energy barriers remain finite while the 
barriers per cross-bridge vanish. The fine structure of the ensuing energy landscape (per cross-bridge) is illustrated 
in Fig. 11(b) for t = t* and successive transitions with a = l/N. At finite N, we see again the two-scale structure, 
however, the macroscopic barrier is markedly higher in a soft than in a hard device. 

In Fig. 12 we illustrate the dependence of the energy landscape in a hard device on the parameter A/, characterizing 
the backbone elasticity in a system with N = 20. To simplify the comparison we adjusted the parameter z = z* at each 
value of A/, so that the coherent states with p — 0 and p = 1 have the same energy (arbitrarily set to zero for the sake 
of comparison). As we know, the snap-spring model converges to the HS model as A/, —» oo and we see in Fig. 12 that 
the “macroscopic barrier” disappears in this limit. The microscopic barriers remain and that is what distinguishes the 
model proposed in (Marcucci and Truskinovsky, 2010) from the HS model. Of course, in the continuum limit N —» oo 
we loose the microscopic barriers and recover the interpolated version of the HS model, see Fig. 5(a). 

In the same figure we illustrate the limit A/, —» 0 and z* —> oo. In this case we recover in the hard device setting the 
predictions of the soft device model with t —» U. Once again, in the continuum limit N —> oo the microscopic barriers 
disappear and we recover the basic picture predicted by the HS model. 

Similar analysis of the barrier structure can be performed for a generic loading. Consider first the hard device 
case. In the continuum limit N —> oo, the variable p becomes continuous. We denote by p,(z) the position of energy 
barrier (saddle point) which is obtained by solving dv/dp = 0 for a given loading z. We obtain 


P*(z) = • 


\[p(A o + A b ) + do A b z - AiAb(z + 1) - Tq/Ii 


y/p(Ao - /li) 

di[di + Ab( 1 + 2z)] - 2vo(/li + Ah) 

2df 


if Ai + Aq, 
if A\ = A 0 . 


Since 0 < p t < 1, the macroscopic energy barrier exists only if < z < Z+, where 



[dz, (do - d0] 1 |^d! - (do + d/,)] 

if d[ * d 0 

Z- = ■ 

(2vo - dj) (/Ij + A b ) 

O 

II 

J4H 


2A\A b 


[A b (do - di)]^ 1 |di (d 0 + A b ) - \fp (do 

+ A b )] if di ± do 

Z+ = ' 

2vo (di + Ab) + A\ (di — Ab) 

if di =d 0 


2AiA b 


Notice that for an arbitrary z, the reaction path does not have to pass through all intermediate configurations, see 
Eq. 18. By inverting Eq. 18, we obtain that for a given elongation, the available metastable states are characterized by 
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Pmf(z) <p< p SU p(z) With, 


p sav (z) — 


Pmf(z) = 


Abz(l - To) - KA o + Ab) 
l(A\ - /lo) + (1 - T 0 )T, 

/lfoz(l - /li) - (l + AO(Ao + Ab ) 

(/ + /li)(/li - /Iq) + (1 - /ll)/ll 


(34) 


Therefore, during the transition the system explores all metastable states satisfying Eq. 34 before cooperatively switch¬ 
ing directly to one of the homogeneous state with p = 1 or p = 0. 



Figure 13: Overall behavior of the snap-spring model. The main frame shows the tension-elongation relations corresponding to p = 0,0.1,..., 1, 
for a system with N = 10 (gray lines) and the tension-elongation relations in the global minimum (thick lines). The dotted lines show the limit 
of the domain where a macroscopic energy barrier is present in a hard device (vertical lines) and in a soft device (horizontal lines). The satellite 
frames (above, hard device; right, soft device) show the energy barriers corresponding to various transitions shown by dashed lines in the main 
frame. In particular, the transition C—>D (hard device) and the transition R—»S (soft device) correspond to z = Z* and t = t*, respectively. The 
other hard device transitions, A-^B and E— >F are outside the bistable interval [z~, z+] and thus contain a collective transition. The other soft device 
transitions, P—>Q and T—>U exhibit the same kind of synchronization along the segment of the ’reaction path’ which is located inside the bistable 
domain. Parameters other than N are as in Fig. 9. 


These results are summarized in Fig. 13 where the main frame shows the tension-elongation relations correspond¬ 
ing to the different metastable states for a system with N = 10. The satellite frames show the transition path related 
landscape corresponding to three different values of z. For instance, in the case where z < Z* (see A—>B) all metastable 
states with p < 7/10 are available and then the system coherently switches from the state with p = 7/10 to the state 
with p — 1. On the main frame, we indicate by vertical dotted lines the values of ’ and z+ which delimit the domain 
where the macroscopic barrier exists. 

A similar analysis can be conducted in the soft device case. First, considering the continuum limit N — > °o, we 
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Dimensional Non-dimensional 


a 

10nm 



i 

-0.8 nm 

/ 

-0.08 

k 

2.7 pN.nnT 1 

N 

100 

kb 

135pN.nm -1 

Ab 

0.5 

k\ 

0.3 pN.nnT 1 

di 

0.23 

h 

k b e 

0.8 pN.nm~' 

3.8 x 10 21 J 

do 

0.44 


Table 1: Realistic parameters for the snap-spring model applied to skeletal muscles. 


solve dw/dp = 0 to locate the position of the energy barrier . We obtain. 


P*(t) 


[(do - d|) -y^lj [(/!] - d 0 ) t + do (d] + -y^/)j 
t~VQ 1 

.d x + 2 


if d, 4 d 0 
if dj = do- 


Then the condition 0 < p,(t) < 1 gives the interval [f , f + ] where the macroscopic barrier is present. Here 


[d() — di] 1 [do(di + -y^/)] 

if do * d 

VO _ di/2 

II 

O 

M-h 

[do - di]^ 1 [d! [d 0 + -y^/)] 

if d 0 * A 

VO + di/2 

II 

O 


Notice that we recover the same boundaries as in the HS model if we consider the limit of infinite stiffness Aj —> 1 
and set di = do. Finally, for a given 1, the available metastable states satisfy p m ((t) < p < p sup (t). The boundaries can 
be obtained by inverting Eq. 23 


Psup(t) — 


t(A o - 1) + /d 0 
/do + (do - / - l)di 


_ f(di - 1) + do(/ + di) 
Pmd j “ /do + (do - / - 1 )d! ' 


Again this shows that (i) the interval where the macroscopic barrier is present is finite and that (ii) the transition 
between the globally stable homogenous states may involve coherent switches depending on the value of the applied 
tension. We illustrate these results in Fig. 13, where the satellite plots on the right show the energy barrier as function 
of £, and where the horizontal dotted lines delimit the macroscopic barrier domain. 

We finally remark that while the individual local minima of the energy landscape at finite N can be interpreted as 
distinct chemical states, the exponentially growth of the number of such states in the thermodynamic limit makes this 
interpretation unpractical. Instead, the description in terms of macro-wells and the associated mechanisms can lead to 
a meaningful quasi-chemical representation of the loading-induced unfolding process. In the case of multiscale un¬ 
folding processes, the macro-wells corresponding to synchronized configurations are usually hierarchically structured 
which makes the mechanical description of such systems an interesting challenge. 


4. Application to skeletal muscles 

In this section, we briefly discuss how our results can be applied to skeletal muscles. Our Table 1 summarizes the 
realistic values of parameters calibrated in (Caruel et ah, 2013). We chose N = 100 by taking into consideration that 
only one third of cross-bridges are attached at any given moment of time (Piazzesi et ah, 2007). 

In Fig. 14 we show the metastability domain (gray area) and the tension-elongation relations describing the global 
minimum of the energy. In the satellites frames we illustrate two sections of the energy landscape along the paths 
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connecting the two homogeneous states: in a hard device at z = z*, see Fig. 14 (a), and in a soft device at t - t,, see 
Fig. 14 (b). These plots are obtained by computing the energy Nv, see Eq. 26 and the energy Nw, see Eq. 31 for hard 
and soft device, respectively. In addition we show the size of the individual energy barriers for single folding (£_,) 
and unfolding ( B,_) events, see Fig. 14 (c and d). The results for finite N (open symbols) were obtained by using 
Eq. 28 for the hard device case and Eq. 32 for the soft device case. In the thermodynamic limit (solid lines) we used 
Eq. 29 in the hard device case and Eq. 33 in the soft device case. Observe that the zone of bistable behavior in a soft 
device, is much broader that in a hard device and spans almost the entire metastability domain. This correlates with 
the fact that the macroscopic energy barrier in a soft device, see Fig. 14(b) is about three times higher than in a hard 
device, see Fig. 14(a). 

Since experiments on muscle fibers are performed at finite temperature, it is of interest to compare the height 
of the microscopic and macroscopic energy barriers with the typical energy of thermal fluctuations kt,6 where kb is 
the Boltzmann constant and 6 is the absolute temperature (see similar analysis of microscopic barriers for titin in 
Benichou and Givli (201 1)). 




B 


B 




z (in nm per half-sarcomere) 


Soft device 



Macroscopic barrier domain in a hard device 


Figure 14: Result of the soft spin model with parameters adjusted to fit experimental data, see Caruel et al. (2013). Main frame, tension-elontation 
relation for a single half-sarcomere. The gray area shows the existence domain of the metastable states. Solid lines, global minimum tension- 
elongation relation in hard and soft devices, (a) and (b), Energy landscape corresponding to single transition between the homogeneous states in a 
hard device (A^B), see (a) and in a soft device (C—>F), see (b). (c) and (d) size of the energy barriers corresponding to the individual folding (B_») 
and (£<_) transitions obtained with both finite N (open symbols) and in the thermodynamic limit (solid lines). To improve readability of the finite 
N plots, we only show one out of five micro-transition points. 


In a hard device, the full range of metastable states between the two homogeneous states (labeled by A and B) 
is available, see Fig. 14 (a). In Fig. 14 (c) we show the size of the microscopic barriers corresponding to single 
folding (resp. unfolding) events, see Il > (resp. B,k). These energy barriers are of the order of ki,f) which suggests that 
thermal noise alone can unfold single bistable elements. However, at small values of £ the height of the energy barriers 
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corresponding to folding is systematically higher than for unfolding and vice-versa for large values of A Therefore, 
although individual folding/unfolding transitions can occur due to thermal fluctuations, the presence of a macroscopic 
barrier, which is at least 50 times higher, ensures that the system is maintained in globally synchronized states, here 
corresponding to points A or B. 

In a soft device, see Fig. 14 (b and d), the situation is slightly different. One can see that at t — t,, the reaction path 
passes through a reduced fraction of metastable states with high energy in the interval 0.3 < £ < 0.8. Therefore, to 
switch from one homogeneous configurations to another (say, from C to F), the system must first reach the metastable 
states D or E through a massive collective transition (C—»D or F—>E) requiring an energy of about 100 k/,6, see Fig. 14 
(b and d). These massive transitions are very unlikely, moreover the reverse transitions have almost zero energy barrier. 
Once the system has reached the metastable states, D or E, we are back to the pattern encountered in the case of a hard 
device where energy barriers for individual transitions were of the order of kb6. The fact that the macroscopic barrier 
in the case of soft device is much higher than in the case of hard device may be the reason behind the anomalously 
slow kinetics of relaxation observed in isotonic experiments (Piazzesi et al., 2002; Reconditi et ah, 2004; Decostre 
et ah, 2005). 

To summarize, individual contractile units in skeletal muscles can be viewed as an assembly of nanometer sized bi¬ 
stable mechanisms. In the absence of long-range coupling, these mechanisms would transform individually producing 
strongly inhomogeneous temporal and spatial microstructures. The presence of long-range interactions is a way for 
the system to strongly bias homogeneous states and in this way passively synchronize individual bistable mechanisms. 
The transition rate associated with a macroscopic transition is extremely sensitive to the number of switching elements 
as the (macroscopic) energy barrier is proportional to N. This suggests that the system can fine-tune its kinetics by 
recruiting a particular number of crosslinkers. Our computation of the individual barrier heights, using the asymptotic 
formulas given by Eqs. 29 and 33, shows that the energy per crossbridge in the asymptotic regime of infinitely large 
N is already well approximated at the realistic values N ~ 100. This means that this number of cross-bridges, 
characterizing a single half-sarcomere, is sufficient for the system to achieve the maximum synchronization ability. 
Even though our conclusions are reached based on the analysis of the internal rather than the free energy, they are fully 
confirmed by the detailed study of the finite temperature effects on muscle dynamics based on the direct modeling of 
the Langevin dynamics (Caruel et al., 2013). 

5. Conclusions 

A parallel bundle of bistable snap-springs is a simple mechanical system which, however, teaches important 
lessons. The presence in this system of mean-field type interactions, induced by the coupling of individual bistable 
units through a common backbone, produces a peculiar mechanical behavior. For instance, the transition between the 
two states (folded and unfolded) takes the form of a collective switching event rather than a sequence of transitions in 
individual elements. The system exhibits negative stiffness and, even in the continuum limit, the mechanical behavior 
is different in soft and hard loading devices. Such systems, where each element is linked with almost equal strength 
with all other elements and the whole is not a sum of the parts are ubiquitous in biology with skeletal muscles 
providing just one of the many examples. While in such systems the individual units may be submitted to random 
thermal fluctuations, the domineering long-range interactions provide a highly efficient way to maintain the individual 
units passively synchronized. More specifically, long-range interactions impose a strong bias between forward and 
backward reaction rates for microscopic transitions favoring globally ordered states. 

The prototypical model studied in this paper has important applications outside the skeletal muscle context. We 
have already mentioned the phenomena of synchronized unzipping of adhesive clusters (Erdmann and Schwarz, 2007; 
Chen and Gao, 2011; Yao and Gao, 2006; Gao et al., 2011; Erdmann et al., 2013) and the cooperative flip-flopping 
of macro-molecular hairpins (Liphardt et al., 2001; Prados et al., 2012; Bosaeus et al., 2012; Woodside et al., 2008). 
A more complex but related example is provided by the “fracture” avalanches during unfolding of macro-molecules 
(Srivastava and Granek, 2013). The broad applicability of the proposed mechanical perspective is also corroborated 
by the fact that proteins and nucleic acids behave differently in isometric and isotonic conditions and that these 
mechanical systems can exhibit negative stiffness (Gerland et al., 2003; Bornschlogl and Rief, 2006; Thomas and 
Imafuku, 2012). In the same vein, the importance of the topology of interconnections among the bonds and the link 
between the cooperativity of unfolding and the dominance of parallel bonding have been long stressed in the studies 
of protein folding (Dietz and Rief, 2008). By emphasizing the crucial role of the force transmitting backbones our 
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study provides a simple paradigmatic description of this class of phenomena. On a more practical side, the model 
suggests an explicit path towards designing bio-mimetic materials and molecular nano-machines whose functioning 
depends essentially on long-range feedback between multi-stable units (Yurke et al., 2000). 
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Appendix A. Stability in a hard device 


The analysis of stability is similar in the HS and the snap-spring models. A technical complication is that in both 
cases the bistable potentials m hs and u ss are singular. Thus, in the HS model the energy wells are infinitely narrow 
and the energy barrier is formally infinite. In the snap-spring model, the energy wells have a finite curvature however 
the spinodal region is reduced to a single point. To study stability, we first remove these singularities by considering 
smoother potentials and then perform the appropriate limiting transition to singular potentials. 

We start with the more general snap-spring model and regularize the bistable potential u ss by introducing an 
extended spinodal interval [Z - e; Z + e] where the new potential S ss is concave. Assume that outside this interval the 
potential u ss coincides with u ss and is therefore convex. In a hard device the new snap-spring energy can be written in 
the form 


v(x,y;z) = 


i N 

iy 

N 


i=i 


S.W+jO’-x,) 


+ f(z-y) 2 . 


(A.l) 


We now analyze stability of the system described by energy (A.l). At a given z, the equilibrium equations can be 
written as 


dv_ 

dxi 


= Mss (*/) + (X ~ )0 = 0, 


i = 1 


(A.2) 


where y(z, p, q, r ) is the equilibrium value of y for a given configuration ( p, q, r ) given by Eq. 14. Assume that within 
the spinodal region there is an interval where w"(x) < -1. Then each of the equations (A.2) has up to 3 solutions. 
We denote these solution by x\ (z, p, q, r), xo(z, p, q, r) and x*(z, p, q. r). The first two solutions x\ and xb correspond 
to the two convex wells of the potential S ss (x) so that i) < l - e and xq > l + e. The third solution x* describes the 
crosslinker in the spinodal region and therefore Z - e < x* < Z + e for all e > 0. We can now compute 


d 2 v(x,y\z) 


dxf 

d 2 v(x,y,z ) 


dx 2 


p,q,r,Xi=xi ,y=y,x&=xj 


p,q,r,Xi=x*,y=y,Xj±i=Xj 


d 2 v(x,y\z) 


<9x7 


p,q,r,Xi=x 0 ,y=y,Xj t i=Xj 


k\ + 1 = hi > 0, 


5"(x«) + 1 


h*(z,p,q, r) < 0, 


ko + 1 = h{) > 0. 


Here h, is negative because the corresponding crosslinker is the spinodal state. The other second derivatives of the 
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energy can be computed explicitly 


d 2 \ix,y,z) 


dxjdxj 

d 2 v(x,y;z) 


dxjdy 
d 2 v(x,y,z ) 


dy 2 


= 0 for i + j, 

p,q,r,x l =x i ,xj=x j ,x ktiy j=x k ,y=y 

= -1, i = 

p,q,r,Xi=Xi,Xj±t=Xj,y=$ 

— 1 + dfo. 

p,q,r,Xi=Xi,y=y 


(A.3) 


To write the expression for the Hessian matrix H (z,p,q,r), it is convenient to introduce the following auxiliary 
quantities 

d 2 \ ix,y,z) 


Hi(p,q,r;z ) = 


dx 2 


p,q,r,x i =x i ,y=y,Xj H =x j 

Each of the variables //, can take three values: h\, ho and h*. Now we can write 


H (p,q,r,z) = 


Hi 0 
0 

0 

v-1 ... 


0 


-1 


■■ 0 : 

0 H n -1 

... —1 1 + Abj 


(A.4) 


To obtain a similar Hessian matrix for the HS model, we need to perform the limit k\o(z) —> °° which means h t —> oo, 

ho — > oo and h, —> -oo, and to drop in (A.4) the last line and the last column. 

From the form of the matrix H, one can see that as soon as one of the terms //, is equal to /;*(’), which means that 
q > 0, at least one of the principal minors of H becomes negative. This means that the absence of crosslinkers in the 
spinodal region is mandatory for stability at e > 0. 

We can now consider the limit e —» 0. The value of the equilibrium strain in the spinodal region x, remains 

between / - e and l + e and thus converges to l, when e —> 0. Therefore, iff/ 0, the configuration ( p,q,r ) is 

necessary unstable and we know that such configurations are necessarily singular. Hence, in our snap-spring model, 
among the (N + 1 )(N + 2) /2 equilibrium branches, N (N + 1) /2 singular branches are unstable which leaves N + 1 
nonsingular branches describing local minima of the energy. Since in the HS model the spinodal states are absent, all 
configurations are automatically metastable. 


Appendix B. Stability in a soft device 

In the soft device case, the tension t is fixed while z becomes an additional degree of freedom. Then the energy of 
the nonsingular snap-spring system reads 



(B.l) 


The analysis of the equilibrium states remains the same and we can similarly define the diagonal terms of the 
Hessian matrix for the energy w{x, y, z\ f). 


Hi(p, q, n t) 


d 2 w(x,y,z, t) 


dx 2 


p,q,r,Xi=Xi,y=y.z=z,Xj*i=Xj 


Here each term can take the following three values h\ — k\ + 1 >0, ho(p, q. r, t) — ko + 1 >0 or li*(p, q, r, t) < 0. The 
other entries are the same as in the hard device case, see Eq. A.3 except that now we have one additional row and one 
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additional column. 


d 2 w(x,y,z\ t) 


dxjdz - - ■ 

p,q,r,x i =5t i ,Xj* i =Stj,y=y,z=p,q,r;t 

for i = 1,..., N 


= 0, 


d 2 w(x, y, z; t) 


dydz 

d 2 w(x,y,z\ t) 


- -Ab, 


p,q,r,Xj=Xj, Xj#=Xj,y=y,z=z 


dz 2 


Ab- 


p,q,r,Xi=Xi,Xft t =Xj,y=yj=z 

By bringing all these second derivatives together we can write the Hessian matrix for the snap-spring model in a soft 
device 

(H l 0 ■■■ 0 -1 0 \ 


H (p, q , r; t) = 


0 


■ 0 : : 

0 ■■■ 0 H n -1 0 

-1 -1 (l+d h ) -A h 

0 0 —Ab Ab y 

A straightforward adaptation of the above analysis shows that, as in a hard device, the system in a soft device is 
unstable only when q + 0, i.e. when at least one cross-bridge is in the spinodal state. 

Finally, to obtain the Hessian matrix for the HS system we need to drop the last row and the last column and 
consider the limit &i,o(z) —» °° which means h\ —> oo, h {) —> oo and h t —> -oo. We also require that A/, —> 0. Then, the 

Hessian reads 

(H\ 0 0 -l\ 


H HS (P, q, r; t) 


0 


-1 


0 : 
0 H n -1 
... -1 1 


The analysis here is similar to the case of a hard device and the conclusion is that again all equilibrium configurations 
are metastable. 

Appendix C. Energy barriers in a hard device 

The explicit expressions for the coefficients in the the asymptotic development (29) are 


a(p',z) = N * 


b(p\z) = N; 


a'(p\z) = N„ 


b'(p\z) = N; 


u ss (l) ~ i'o + 


\\A h + Mp) 

/ - (A b z - pdi)} 

-d 0 [A b z- pA\} 2 

2 

Ab + A(p) 

r 


{[ d/, + Mp)] / - (1 - A 0 )(A b z - pd|)} 2 


2 |^/1& + A(p)J 


Mss(0 + 


[[Ab + A(f)] (/+D — [A h (z + 1) + (1 - p)A 0 ] 

} di [Ab(z + 1) + (1 - p)d 0 ] 2 

2 

Ab + A(p)]‘ 



{[d* + M P )\ (l + 1) - (1 - dj) [A b (z + 1) + (1 - ,p)d 0 ]} 2 


2 [/tfc + A(/?)] 


with A (p) given by Eq. 25. 
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Appendix D. Energy barrier in a soft device 


The explicit expressions for the coefficients in the the asymptotic development (33) are 

(A (p)Z-a-M)' 2 

«ss(0 - >'o + 


c(p; t ) = N, 

d(p\ t ) = N 2 
c'(p ; t) = N* 

d'{p-t) = Nl 


(Mp)l ~ (t - pAi)) - A()(t - pA\) 2 


2A (p) 2 


(Mp)I - (1 - d 0 )(f - pAi)) 2 


UssO) + 


2k(p)i 

(a mi + 1) - (t + (1 - p)A 0 )f - A 1 [t + (1 - p)A 0 ] 2 
2A(P) 2 


(A(p)(Z + 1) - (1 - Ai)(t + (1 - p)A 0 )f 
2A(P) 3 


with A (p) given by Eq. 25. 
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